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ABSTRACT 

Context. A molecular survey recently performed ahead of HH 2 supports the idea that the observed molecular enhancement is 
due to UV radiation from the HH object. 

Aims. The aim of the present work is to determine whether all HH objects with enhanced HCO + emission ahead of them also 
exhibit the same enhanced chemistry as HH 2. We thus observed several molecular lines at several positions ahead of five 
Herbig-Haro objects where enhanced HCO + emission was previously observed. 

Methods. We mapped the five Herbig-Haro objects using the IRAM-30 m. For each position we searched for more than one 
molecular species, and where possible for more than one transition per species. We then estimated the averaged beam column 
densities for all species observed and also performed LVG analyses to constrain the physical properties of the gas. 

Results. The chemically richest quiescent gas is found ahead of the HH 7-11 complex, in particular at the HH 7-11 A position. 
In some regions we also detected a high velocity gas component. We find that the gas densities are always higher than those 
typical of a molecular cloud while the derived temperatures are always quite low, ranging from 10 to 25 K. The emission of 
most species seems to be enhanced with respect to that of a typical dense clump, probably due to the exposure to a high 
UV radiation from the HH objects. Chemical differentiation among the positions is also observed. We attempt a very simple 
chemical analysis to explain such differentiation. 


Key words. ISM: Herbig-Haro objects — ISM: abundances — ISM: clouds — ISM: molecules — Radio lines: ISM — Stars: 
formation 


1. Introduction 

Herbig-Haro (HH) objects are regions where protostel- 
lar jets interact hydro-dynamically with the surround¬ 
ing molecular cloud, resulting in strong atomic line emis¬ 
sion characteristic of shocks. There is now ample and in¬ 
creasing evidence that HH objects also have a less direct 
but very diverse chemical influence on nearby molecu¬ 
lar clumps which have not yet been reached completely 
by the flow itself. In fact, the radiation generated in 
the HH objects, as well as the shocks produced by the 
HH objects alter the chemical composition of the molec¬ 
ular gas surrounding them in a very complex way. An 
example of the complex interaction of the HH object 
with its surroundings, and the most studied one, is the 
HH 2 region (e.g. IGirart et al. 20021 IDent, et al. 20031 
ILefloch et al. 20051 IGirart et al. 20051) . Some of the gas 
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ahead of the HH object seems to be altered solely by 
the UV radiation: in the last decade or so quiescent, 
cold, molecular condensations have been found ahead 
of several HH objects. These cool condensations were 
found in enhanced emission of HCO + and NH 3 , which 
are tracers of high d ensity gas (e.g. fRudolph et al. 19881 
ITorrelles et al. 19921 IGirart et al. 19941) . These" emissions 
probably arose because the UV radiation generated in 
the HH shock evaporated the icy mantles on dust grains 
in small density enhancements in the cloud near to the 
head of the jet (IGirart et al. 199411 . Theoretical and fur¬ 
ther observational studies (IViti fc Williams 19991 here¬ 
after VW99; IViti et al. 20031 and IGirart et al. 200511 con¬ 
firmed this picture and showed that, because of the quies¬ 
cent nature of such clumps (clear from the linewidths of 
the emission as well as the derived kinetic temperature of 
the gas), the high abundances of HCO + and other species 
could not be due to shock chemistry. 
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In this paper, we are particularly interested in these 
anomalous, quiescent enhanced molecular emissions from 
localized regions in the vicinity of the HH objects. The 
VW99 models succeeded in explaining the known enhance¬ 
ments of HCO + and NH 3 relative to CS, predicted the 
amount by which these species should be increased, but 
also indicated other species which should be enhanced 
(by factors of 10-200 with respect to dark molecular 
clouds e.g L1554), in particular CH 3 OH, H 2 S, SO, S0 2 , 
H 2 CO, C 3 H 4 , H 2 CS. Such predictions were confirmed 
by recent observations and further modelling of HH 2 
HGirart et al. 20021 lYiti et al. 20031 IGirart et al. 200511. It 

is important to note however that while single dish ob¬ 
servations HGirart, et al. 20021 1 Viti et al. 200311 show that 
the chemistry of the illuminated quiescent clumps ahead 
of the HH 2 object is essentially consistent with the UV 
field from the HH 2 front releasing the grain mantles 
and hence enhancing the gas chemistry, BIMA maps of 
such region in several molecules expected to be enhanced 
IlGirart et al. 20051 have shown how the HH 2 object af¬ 
fects the chemistry and dynamic of the surrounding envi¬ 
ronment in different degrees. In summary, although Girart 
et al. (2005) concluded that a very complex morphological 
kinematical and chemical structure of the molecular gas is 
present ahead of the HH 2 object, they also confirmed that 
there is indeed an apparently quiescent molecular compo¬ 
nent that is more exposed to the UV radiation from the 
HH 2 front; the analysis of observational molecular lines 
from this region can be explained by UV radiation from 
the HH shock lifting the ice mantles from grains of a multi¬ 
density component clump and the consequent chemistry 
that arises is as complex and varied as indicated in the 
VW99 models. 

The question we attempt to answer in this paper is the 
following: is the rich photochemistry characteristic of UV 
illumination peculiar to the HH 2 or is it common to all 
HH objects? The simplest and most direct way to answer 
such question is to perform a survey of molecular species 
predicted to be enhanced by the VW99, ahead of several 
HH objects where enhanced quiescent HCO + and/or NH 3 
have already been detected. We have now performed such 
a survey with the IRAM 30-meter and we report the ob¬ 
servational results in this paper. The main aims of this 
paper are: 1) determine whether other HH objects exhibit 
the same peculiar enhanced chemistry of HH 2 as pre¬ 
dicted by VW99 and 2) assuming the HH object to be the 
source of the UV radiation, look for correlations between 
the enhanced molecular emission and the distance of such 
emission to the HH object. 

We report our observations in § 2; in § 3 we present 
our results and the derivation of beam averaged column 
densities. In § 4 we analyse our data by means of radiative 
transfer modelling. In § 5 we briefly discuss the chemical 
characteristics of the gas and in § 6 we summarize our 
findings. 


2. Observations 


The observations were carried out with the IRAM 30 m 
radio telescope on Pico Veleta in November 2000. Our 
sample consisted of five Herbig Haro regions; for most of 
them we took several positions where enhanced HCO + 
was previously detected (see individual subsections and 
Tabled. For each position we searched for several molecu¬ 
lar emissions; when possible we searched for more than one 
transition of the same species (see Table 0 . The molec¬ 
ular species chosen are those which were expected from 
VW99 models to be enhanced (with respect to dark molec¬ 
ular clouds) by reasonably large factors, along with trac¬ 
ers of the underlying column density in the cloud (CS 
and CO isotopes) and tracers of the ionization (the ra¬ 
tio DCO + /HCO + : ICaselli et al. 19981) . Due to time con¬ 
straints, we could not attempt to detect all the species 
predicted by the VW99 at every position. Our criteria 
varied from object to object but in general if HCO + emis¬ 
sion was weak we did not attempt a detection in other 
species. We also prioritized species easier to observe, so 
for example, although both SO and S0 2 are predicted by 
the models to be enhanced we concentrated on SO as its 
lines are usually stronger and also it has a higher pre¬ 
dicted abundance. All the transitions were observed in 
the frequency bands of 100, 150, 230 and 270 GHz (see 
Table 2 for line frequencies). Integration times were typ¬ 
ically 2-12 minutes (position switched), with longer at¬ 
tempts at detecting weak lines (e.g. NS: 30 mins). The 
spectral resolution was 40 kHz (100/150 GHz bands) or 
80 kHz (230/270 GHz bands) resulting in a velocity reso¬ 
lution typically of 0.05 km/s and at most 0.13 krn/s. The 
bandwidth was at least 20 MHz (100/150 GHz bands) 
or 40 MHz (230/270 GHz bands). Conditions were gen¬ 
erally poor and variable for winter nights with l- 6 mm 
PWV (precipitable water vapour). Resulting system tem¬ 
peratures were Tsys ~ 100-120, 200-460, 200-300 and 
400-700 K for the 100, 150, 230 and 270 GHz bands re¬ 
spectively, depending on frequency and weather. Pointing 
was good to within 4 

arcsec. Spectra were corrected for the beam efficiencies 
B eS = 0.75, 0.69, 0.51, and 0.42 in the 100, 150, 230, 
and 270 GHz bands respectively. The 30m beam FWHM 
ranges as i / _1 from 28" at 86 GHz down to 11" at 230 GHz. 
See on-line Tables El to E3 for details of the molecular 
transitions observed. 


The choice of our sample of objects was made by se¬ 
lecting all the HH objects with known HCO + emission 
enhancements. For some objects, we detected two velocity 
components in several species: a high velocity one, proba¬ 
bly coming from the closest part to the HH object where 
the gas has already been dynamically affected; and an am¬ 
bient one. In this paper we will mainly concentrate on the 
latter. 
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Table 2. Frequency setups of the IRAM 30m observations and list of positions observed. 


Molecular 

Transitions 

Frequency 

(GHz) 

Beam 

(") 

HH 7 

-11 



HH 34 


N2264G a 

HH 1 

A 

X 

B 

A 

E 

c 

D 

1 

2 

13 CO 2-1 

220.3980 

11.0 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

C ls O 2-1 

219.5604 

11.0 

X 

X 


X 

X 



X 

X 

X 

HCO+ 1-0 

89.1885 

28.0 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

HCO+ 3-2 

267.5576 

9.0 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

H 13 CO + 1-0 

86.7543 

29.0 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

H 13 CO+ 3-2 

260.2555 

9.5 

X 

X 

X 

X 

X 



X 

X 

X 

DCO+ 2-1 

144.0773 

17.0 

X 

X 


X 

X 



X 

X 

X 

DCO+ 3-2 

216.1126 

12.0 

X 

X 


X 

X 



X 

X 

X 

HCN 3-2 

265.8864 

9.0 

X 

X 


X 

X 



X 

X 

X 

H 2 CO 2 0 ,2-1o,i 

145.6029 

17.0 

X 

X 


X 

X 





X 

H 2 CO 31,3-2 !, 2 

211.2115 

12.0 

X 




X 






CH 3 OH 2 k -l k b 

96.7414 

26.0 

X 

X 


X 

X 



X 

X 


CH 3 OH 3fc-2fc c 

145.1032 

17.0 

X 

X 


X 

X 



X 

X 


CH 3 OH 5fc-4fc d 

241.7914 

10.0 

X 

X 

X 








C 3 H 4 5o~4o 

85.4573 

29.0 

X 




X 






C 3 H 2 2i, 2 -1o,i 

85.3389 

29.0 

X 




X 






CsH 2 3i, 2 -2 2 ,i 

145.0896 

29.0 

X 

X 









CsH 2 4i,4-3o,3 

150.8519 

17.0 

X 










CS 3-2 

146.9690 

17.0 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

CS 5-4 

244.9356 

10.0 

X 




X 






SO 3 2 -2i 

99.2999 

20.0 

X 

X 


X 

X 



X 

X 

X 

SO 65-54 

219.9494 

11.5 

X 

X 


X 

X 



X 

X 

X 

SO 7 6 -6 5 

261.8437 

9.5 

X 

X 


X 

X 



X 

X 

X 

H 2 S li,o-0 0 ,i 

168.7628 

14.5 

X 

X 

X 





X 

X 

X 

H 2 CS 4 0 , 4 -3o,3 

137.3711 

18.0 

X 




X 






HCS+ 5-4 

213.3606 

12.0 

X 










NS 7/2—5/2 

161.2972 

15.0 

X 




X 






OCS 7-6 

85.1391 

29.0 

X 











a Abbreviation of NGC 2264G JMG1 and JMG2 
b 2o-lo A + transition. 
c 3o-2o A + transition. 
d 5q-4o A+ transition. 


Table 1. Source list 


Object 

RA 

J2000 

DEC 

J2000 

V LSR 
kms -1 

HH 7-11 A 

03:29:06.5 

31:15:35.5 

+ 8.0 

HH 7-11 B 

03:29:06.9 

31:15:47.4 

+ 8.0 

HH 7-11 X 

03:29:05.7 

31:15:30.5 

+ 8.0 

HH 1 

05:36:20.2 

-06:45:01.3 

+9.0 

HH 34 A 

05:35:34.3 

-06:29:07.0 

+8.5 

HH 34 C 

05:35:32.2 

-06:29:04.9 

+8.5 

HH 34 D 

05:35:32.0 

-06:29:23.9 

+8.5 

HH 34 E 

05:35:30.6 

-06:29:03.8 

+8.5 

JMG 99 Cl 

06:41:04.5 

09:56:02.0 

+4.6 

JMG 99 G2 

06:41:07.1 

09:56:07.0 

+4.0 


2.1. Brief description of the sources 

2.1.1. HH 7-11 

The HH 7-11 system (d ~220 pc: ’ Cernis 1990 ) consists 
of fairly low excitation HH objects, with HH 7 being at 


the far end of the outflow IjSolf et ah 198711 . However, the 
region surrounding HH 7-11 is very rich in young stel¬ 
lar objects (e.g. Rodriguez et al. 1999), hence the to¬ 
tal radiation impinging on the surroundings is unknown. 
Enhanced HCO + 1-0 emission has been previously ob¬ 
served just ahead of HH 8-10 objects in the form of small 
clumps with 7" sizes (1500 AU, |Rudolph et al. 1988| ). CS 
and NH 3 ( [Rudolph et al. 2001| are also reported towards 
these clumps. Submillimeter observations in the region re¬ 
veal a dust source, HH 7 11 MMS 4, that engulfs the 
HCO + clumps ifChini et al. 200 1). From the dust observa¬ 
tions Chini et al. (2001) derive that the total mass of the 
quiescent gas around HH 8-10 is 4.0 M 0 . From the study 
of the HCO + 3-2 emission, Dent et al. (1993) suggest that 
the dense molecular gas and dust around HH 8-11 is the 
densest, slowest component of the outflow. This is some¬ 
what similar to the picture of the molecular gas ahead of 
HH 2 (IGirart et al. 200511 . 
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Table 3. Averaged beam column densities of the low velocity component in the clumps ahead of HH objects observed 
with IRAM 30m (Columns 2-10); the last column shows, for reference, the column densities (when available) at position 
13 ahead of HH2 taken from Girart et al. (2005, table A.2). To fit the table, we abbreviated HH 7-11 to 7-11, the HH 
34 clumps to HH34, and NGC 2264 G JM G2 to JMG 99 G2. a(b) stands for axlO h . 


MOL. 

7-11A 

7-11X 

7-11B 

HH34 A 

HH34 E 

HH34 C 

HH34 D 

HH 1 

JMG 99 G2 

HH 2 

CO“ 

2.6(18) 

2.3(18) 

2.1(18) 

1.9(18) 

1.4(18) 

1.3(18) 

1.2(18) 

3.5(18) 

3.7(17) 

3.3(17) 

^co 

4.6(16) 

4.0(16) 

3.7(16) 

3.0(16) 

2.2(16) 

2.0(16) 

1.9(16) 

5.5(16) 

6.5(15) 


CH 30 H 

3.5(14) 

4.1(14) 


2.9(13) 

4.9(13) 



5.6(13) 

9.8(13) 

6.2(13) 

HCO+" 

2.0(14) 

1.6(14) 

2.2(14) 

1.2(13) 

3.3(13) 

1.7(13) 

9.7(12) 

2.1(13) 

5.0(12) 

4.0(13) 

H la CO+ 

2.3(12) 

1 .8(12) 

2.5(12) 

1.5(12) 

4.2(12) 



2 .6(12) 

6 .0(11) 


SO 

8.5(13) 

8.6(13) 


8 .2(12) 

1.7(13) 



1.7(13) 

1.6(13) 

1.7(13) 

H 2 CO c 

9.3(13) 

9.6(13) 


1.0(13) 

1.1(13) 



1.6(13) 


9.0(13) 

cs 

6.0(13) 

6.2(13) 

5.9(13) 

6 .8(12) 

8 .2(12) 

6 .2(12) 

5.6(12) 

1.0(13) 

1.1(13) 

2 .0(12) 

h 2 s 

5.0(13) 

4.7(13) 

4.5(13) 





1.2(13) 

6 .0(12) 


C 3 H 4 

5.7(13) 




<2(13) 






HCN 

2.3(13) 

1.6(13) 


<6(11) 

<6(11) 



1 .8(12) 

7.8(12) 

2 .8(12) 

OCS 

2.1(13) 










H 2 CS 

1.1(13) 




<2(12) 






C 3 H 0 

8.9(12) 

8.4(12) 



8.5(11) 






NS 

7.8(12) 




<5(12) 






DCO+ 

5.0(12) 

3.6(12) 


3.8(11) 

2.9(11) 



1.9(11) 

2 .2(11) 

7.8(11) 

HCS+ 

<2(12) 











a Optical depth of HH 7-11 B adopted from HH 7-11 X, and of HH 34 C and HH 34 D adopted from HH 34 E. 
6 Optical depth of HH 34 C and HH 34 D adopted from HH 34 E. 
c Adopted an ortho/para ratio of 1.5 



Fig. 1. Sketch of the HH 7 -11 region. The gray scale shows 
the optical R band image of the HH 7 11 objects, while 
the contour plots shows the BIMA HCO + 1-0 blueshifted 
high velocity emission (from Mehringer 1996). The optical 
image was taken from the public catalog Aladin. The cir¬ 
cles are the IRAM 1mm (dotted) and 3mm (solid) beams 
at the A, B and X positions. 


We selected three positions, two of them, 
HH 7-11 A and B, toward the HCO + 1-0 clumps 
I jRudolph et al. 1988| ) and one, HH 7-11 X, toward 


a HCO + 3-2 emission peak southwest of HH 8 
( Dent et al. 19931 . Note that the 3 mm beams of 
our 3 positions (A, X and B) intersect (see Figure^!. 


2.1.2. HH 1 


The HH 1-2 system (d ~450 pc) is one of the most in¬ 
tensely studied HH object system to date. The proper 
motions of HH 1 and HH 2 show that they are moving in 
opposite directions plerbig & Jones 198ip . Their power¬ 
ing source is probably the embedded luminous young star 
known as HH 1 2 VLA 1 (jPravdo et al. 19851) . Infrared 
and optical images show very complex structures, quite 
randomly positioned in the case of HH 2 and more like a 
set of bow shock features in HH 1. Clumpy molecular emis¬ 
sion ahead of HH 2 has been the subject of an extensive 
recent study IlGirart et al. 20051 . Enhanced HCO + and 
NH 3 have also being previously observed ahead of HH 1 
l|Davis et al. 199f)l ITorrelles et al. 19931) . No dust is de¬ 
tected towards the NH 3 clumps (IChini et al. 200 1) up to 
0.22 Jy beam -1 (3-u) at 850 /im. Note that the dust peak 
emission ahead of HH 2 is 0.20±0.04 Jy IIDent et al. 2003|l . 
Our IRAM 30 m observations pointed toward the NH 3 
(1,1) ammonia peak ahead of HH 1 (see Figure El ■ We 
did a search for the most abundant species towards this 
region. 
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Fig. 2. Sketch of the HH 1 region. The gray scale shows 
the optical R band image of the HH 1 object, taken from 
the public catalog Aladin. The circles are the IRAM 1mm 
(dotted) and 3mm (solid) beams. 

2.1.3. HH 34 

The HH 34 is a high excitation large bow shock 
( [Reipurth et al. 2002| ), which is part of a parsec scale HH 
system located in the molecular cloud L1641 in Orion 
(d~480 pc). Its source may be a faint star at 46,000 
AU away. Enhanced, and quiescent, HCO + 1-0 emission 
was detected from five clumps with sizes of about 15" or 
7000 AU |Rudolplr et al. 1992j ). The clumps are located 
south of HH 34 at a distance to the HH object of 0.04-0.09 
pc, corresponding to 8000-18000 AU. We have searched 
for the predicted set of molecules at four positions which 
coincide with the peak HCO + emission of four of the de¬ 
tected clumps, A, E, C and D, by Rudolph et al. (1992; 
see Figure |oJ, although for these two last positions only 
13 CO, HCO + and CS were observed. 

2.1.4. NGC 2264G 

NGC 2264G is associated with a spectacular CO out¬ 
flow driven by a Class 0 protostar { Lada et al. 19961 
IGomez et al. 19941 |Ward-Thompson et al. 19951 ) • Girart 
et al. (2000) found two small clumps (sizes of ^ 14") 
of apparently quiescent HCO + 4-3 emission, JMG 99 
G1 and JMG 99 G2, west of the protostar and just 
ahead from shock-excited near-IR emission. Recent ob¬ 
servations shows that these clumps are also associated 
with strong high velocity SiO and CH 3 OH emission 
I jGaray et al. 20021 . We observed these two clumps. 

3. Results 

Figure O shows examples of the spectra of the molecular 
transitions for which emission was detected, for HH 7 



Fig. 3. Sketch of the HH 34 region. The grey scale shows 
the optical R band image of the HH 34 object, taken from 
the public catalogue Aladin. The circles are the IRAM 
1mm (dotted) and 3mm (solid) beams at the HH 34 A, C, 
D and E positions. 



Fig. 4. JCMT contour maps of the integrated HCO + 4-3 
emission (from Girart et al. 2000) overlapped with the 
grey scale image of the near-IR H 2 emission. The circles 
are the IRAM 1mm (dashed) and 3mm (solid) beams at 
the JMG99 G2 and G1 positions. 

11 A and HH 34 E. HH 7 11 is the region where molec¬ 
ular lines are brighter so position HH 7 11 A was se¬ 
lected to carry out the most complete survey. Apart from 
HCS + , all the molecular species predicted by the VW99 
to be enhanced were detected toward this position. All 
the positions observed but JMG 99 G1 shows that the 
emission is dominated by relatively narrow lines (mostly 
between 0.4 and 1.5 km s -1 ), typical of quiescent emis¬ 
sion, although line width can vary from region to region 
or between different species (see on-line Tables □ □□□ 
and [TUI). HH 1 shows the largest differences between line 
widths for different species, ranging from ~ 0.5 km s -1 
for DCO + to 2.4 km s -1 for CH 3 OH. These large vari¬ 
ations are also observed in the HH 2 molecular clumps 
HGirart, et al. 20021) . For HH 7-11 positions and JMG 99 
G2 most line widths are in the 0.8 to 1.4 km s -1 range. 
Interestingly, for these positions the HCN 3-2 line emis- 
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sion is significantly broader than the averaged line width 
(roughly by a factor 2). 13 CO lines are in most positions 
either clearly non-Gaussian or have several components, 
which can be explained from contribution from lower den¬ 
sity components along the line of sight. In these cases, in 
order to estimate the integrated line intensity, we selected 
the velocity range where the higher density tracer showed 
up. Strong absorption is observed at the center of the 
HCO + 1-0 line in the HH 7-11 and specially in the NGC 
2264G region, which is due to the presence of a low density 
molecular gas in the foreground ( iGirart et al. 2000ll . 

Some of the observed positions show high velocity 
wings in the HCO + 1-0 line (the three positions in HH 7 
11, HH 1, HH 34 E and JMG 99 G2). This is more clearly 
seen in HH 7-11 A (see Fig. 0, where other strong lines 
also show wings: 13 CO, H 2 CO, CS, CH 3 OH. 

The properties of the molecular emission toward JMG 
99 G1 is quite different from the other regions. Some 
of the molecular transitions have emission at the sys¬ 
temic velocity of the cloud (e.g. C 18 0 2-1, HCO + 4-3, 
CS 3-2) but the strongest emission arises from a broad, 
Av cs 5 km s , blueshifted component. This compo¬ 
nent is also observed in SiO l |Garay et al. 2002| . All these 
features indicate that most of the molecular emission in 
JMG 99 G1 comes from shocked outflow material (see 
Section 5.3). 

4. Analysis 

The excitation temperature and the optical depth for the 
HCO + can be derived from the line ratio of the different 
HCO + and H 13 CO + observed transitions. From this anal¬ 
ysis we found that T ex —6-8 K and r(HCO + 3-2)~3-ll. 
Therefore HCO + beam averaged column density is derived 
adopting an excitation temperature of 7 K and correcting 
for the optical depth. For species which also have a high 
dipole moment (e.g. HCN, H 2 CO) the column densities 
are estimated assuming LTE and adopting T ex = 7 K but 
also assuming that the emission is optically thin. For the 
other species we adopt T ex = 10 K. Column densities of 
CO and isotopes are corrected for optical depth. For a non 
detection we report the 3-er upper limit. Table[3]shows the 
derived beam averaged column densities for the quiescent 
emission (i.e. excluding the high velocity component) for 
all the positions but JMG 99 Gl. For reference, we also 
include the column densities in the S0 2 clump ahead of 
HH 2 taken from Girart et al. (2005). Note that Girart et 
al. obtained column densities for several positions along 
the slice in the S0 2 clump (see their Table A.2); we chose 
to present only one of them (position 13, T ex = 13± 2 
K). Most of the lines observed toward JMG 99 Gl show 
emission that is more characteristic of strongly perturbed 
gas. In addition we used the LVG and RAD EX codes to 
carry out a more detailed analysis for the molecules with 
several transitions observed (see Table 4). 

For the high velocity emission toward HH 7-11 and 
JMG 1 we also estimated the column densities assuming 
LTE and adopting T ex = 30 K. In these cases, we com¬ 


puted the column density in different velocity intervals 
along the velocity range here the high velocity is detected 
(see §5.3). 

4.1. LVG and RADEX analysis 

In this section we estimate, to a first approximation, 
some of the physical and chemical characteristics of most 
clumps, such as the density and the temperature by means 
of large velocity gradient (LVG) and RADEX modelling. 

We observed multitransition observations of HCO + 
and H 13 CO+ for HH 7-11 A, HH 7-11 B and HH 7 11 X, 
HH 1, HH 34 A and HH 34 E. We therefore carried out ra¬ 
diative transfer calculations using the LVG model approx¬ 
imation. The main goal of the modelling is to constrain 
the range of densities and temperatures of the clumps, 
as well as the HCO + true column densities. LVG mod¬ 
elling is performed varying three parameters: the column 
density of HCO + , the temperature, and the volume den¬ 
sity of the gas; one parameter is provided as an input, 
and the aim is to obtain a set of LVG solutions in the 
other two parameters plane. For HCO + , we ran the LVG 
code with the HCO + column density as input parame¬ 
ter and we stopped when a set of reasonable solutions 
was obtained. To constrain the set of solutions we used 
all the observed line intensity as well as the observed 
line temperature ratios ([H 12 CO + 3-2]/[H 13 CO + 3-2] and 
[H 13 CO + 2-l]/[H 13 CO + 3-2]). To account for the differ¬ 
ence in beam sizes between the different transitions, we 
multiplied our intensities and temperature ratios by the 
appropriate filling factor (assuming that the size of the 
emission is the same for all the transitions and that the 
emission has a Gaussian shape). We initially used the 
HCO + size derived from previous observations (see fol¬ 
lowing subsections), although we also explored different 
values around the reported values to check if the LVG fits 
improved. The solutions were searched in the Tk~n( H 2 ) 
plane for different HCO + column densities ranging from 
lxlO 13 to lxlO 15 cm -2 (the HCO + beam averaged col¬ 
umn densities for the clumps where we perform the LVG 
analysis are in the 2xl0 13 -2xl0 14 cm -2 range). The re¬ 
sults are summarized in Table 01 The range of values of 
the derived densities and temperatures from the HCO + 
LVG analysis was used to perform the same analysis with 
CS for HH 7 11 and HH 34 E, using as completely free 
parameters only the CS column density and the size of its 
emission. 

SO and C 3 H 2 are the only other species for which, for 
some objects, we observed multiple transitions and colli- 
sional rates are available. For the SO (HH 7 11 A and 
HH 7 11 B) and C 3 H 2 (HH 7-11 A) observed transitions 
we used RADEX, a non-LTE molecular radiative transfer 
code which assumes an isothermal homogeneous medium 
(ISchoier et al. 20051 . To find the possible set of solutions 
for the SO and C 3 H 2 we used the three observed transi¬ 
tions as well as their line ratios with the line intensities 
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Fig. 5. Spectra of selected transitions of the different species observed toward (from left to right) HH 7-11 A, HH 7 
11X, HH 1, HH 34E, HH 34A, JMG 1 and JMG 2. 13 CO 2- 1 spectra is scaled down by a factor of 3 to fit in the boxes. 
The transition lines showed are (from left to right and from top to bottom): 13 CO 2-1, CS 3-2, HCO + 1-0, HCO + 
3-2, SO 3 2 -2i, DCO+ 2-1, H 2 CO 2 0 , 2 -1 0 > i, HCN 3-2, CH 3 OH 2^-1^, and H 2 S li i0 -0 0 ,i. The dashed lines in the CO 
and HCO + panels are for the C 18 0 and HC 13 CO + isotopes at the same transition, respectively. 


Table 4. Best fits with LVG and RAD EX for column densities, gas densities and temperatures. a(b) stands for axlO 6 . 


HCO+ 

Object 

Size 

V(HCO + ) (cm"") 

n( H 2 ) (cm a ) 

T k ( K) 

HH 7-11 A 

20 " 

5(14)-1(15) 

8(4)-3(5) 

15-24 

HH 7-11 B 

20 " 

5(14)-1(15) 

9(4)-4(5) 

15-24 

HH 7-11 X 

20 " 

5(14)-1(15) 

3(5)-2(6) 

7-11 

HH 1 

10 " 

5(14)-1(15) 

4(5)-l(6) 

8-11 

HH 34 A 

10-15" 

5(14)-1(15) 

3-4(5) 

7-11 

HH 34 E 

10 -12" 

5(14)-1(15) 

1.5-9(5) 

8-12 

CS 

Object 

Size 

JV(CS) (cm-") 

n(H 2 ) (cm a ) 

T k ( K) 

HH 7-11 A 

30" 

5(14)-1(15) 

l(5)-3(5) 

15 

HH 34 E 

30" 

3(13)-1(14) 

3(5)-9(5) 

10-12 

SO 

Object 

Size 

N( SO) (cm"") 

n(H 2 ) ( cm d ) 

T k ( K) 

HH 7-11 A 

20 " 

3(14)-1(15) 

3(4)-4(5) 

>10 

HH 7-11 X 

20 " 

3(14)-1(15) 

2(4)-2(5) 

>10 

Object 

Size 

V(C 3 H 2 ) (cm’") 

n(H 2 ) ( cm d ) 

T k { K) 

HH 7-11 A 

20 " 

3(13)-1(14) 

4(4)-2(5) 

>10 


corrected by the filling factor. The set of best solutions 4.1.1. HH 7—11 
are shown in Table Q] 

For all the clumps we find good fits for HCO + column 
densities ranging from 5xl0 14 to lxlO 15 cm -2 . As an 
example of the LVG modelling, Figure|Cl(top) shows a set 
of LVG solutions in the 7Vn(H 2 ) plane for HH 7- 11 A. 
The best fit is obtained for a gas density of ~ 8xl0 4 - 
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3xl0 5 cm -3 , a T% of 15-24 K and a source size of ~ 20"; 
while for HH 7-11 B we find similar solutions, for HH 7 
11 X we obtain a higher gas density range and definitely 
a lower kinetic temperature. Although HH 7-11 X seems 
close to HH 8 , the low temperature may be an indication 
that this clump is further away from the highest excitation 
HH object. Also note that we do not know the line-of-sight 
component of the distance between clump and HH object. 
The derived source size is larger than the estimated size 
from Rudolph & Welch (1988) maps. This may be due to 
the clumps being very close to each other causing their 
emission to be ‘contaminated’, yielding to an apparently 
bigger size. 

There are two CS transitions detected toward HH 7 
11 A, so the LVG analysis have been performed by us¬ 
ing as a constrain the lower limit of the kinetic temper¬ 
ature derived from the LVG HCO + analysis, 71- =15 K. 
Interestingly, we do not find any solution for a emitting 
size less than 30", suggesting that the CS emission is more 
extended than the HCO + one. Adopting a size of 30", 
the set of possible solutions spread out over a large range 
of volume and column densities. Table 0] shows the CS 
column density assuming a volume density similar to the 
value derived from the HCO + . 

For HH 7-11 A and X we also performed RAD EX anal¬ 
ysis of SO and, for HH 7-11 A, of C 3 H 2 . As for the other 
molecules we find that the best SO column density fits 
are reached for values higher than the beam averaged ones. 
From SO and C 3 H 2 we were not able to constrain the tem¬ 
perature of the gas but we notice that it has to be higher 
than 10 K (see Figurecentral and bottom panel). The 
derived gas densities are slightly lower than those derived 
from the HCO + analysis. 

4.1.2. HH 1 

We did not detect the H 13 CO + 3-2 line, hence we 
had to adopt a 3-<r upper limit for it as well as for 
the [H 12 CO+3 - 2]/[H 13 CO+3 - 2] and [H 13 CO+1 - 
2]/[H 13 CO + 3 — 2] line ratios. We find that, unlike for 
HH 7-11 clumps, the best fits are obtained for a source 
size of 10 " consistent with the size found in the literature 
IITorrelles et al . 1 994ft . This may support the idea that the 
apparent larger size of the HH 7-11 clumps is merely a 
consequence of contamination. While the HCO + column 
and gas density of the HH 1 clump are similar to those for 
HH 7 11, its temperature is certainly lower and somewhat 
better constrained. 

4.1.3. HH 34 

We find a remarkably small range of solutions for the ki¬ 
netic temperature, at 7-11 K for HH 34A and 8-11 K for 
HH 34 E. As for HH 1, we find that the best solutions are 
obtained for 10"—15" sizes, consistent with those found in 
the literature | |Rudolph et al. 1992| ). The HCO + column 
densities range from 5xl0 14 to lxlO 15 cm -2 . For HH 



_l_1_i_ 

10 20 30 


T ki „ (K) 

Fig. 6. Plot of the set of LVG and RADEX solutions 
in the Tk-u(J12) plane for HH 7 11 A assuming a source 
size of 20" for the following molecules: Top panel HCO + , 
derived using LVG for a column density of 10 15 cm -2 . The 
grey scale and the thick solid line in its center shows the 
HCO + 3-2 to H 13 CO + 3-2 line ratio. The three solid lines 
show the H 13 CO + 1-0 to H 13 CO + 3-2 line ratio. The thick 
solid and dashed grey line show the the HCO + 3-2 and 
HCO + 3-2 lines, respectively. Central panel SO, derived 
using RADEX for a column density of 3 x 10 14 cm -2 . The 
thin contours levels are (from lowest to highest contour) 
X 2 = 20,16,12, 8 ,4, 2,1. The thick grey contour shows the 
68 % confidence region of y 2 . Bottom panel. C 3 H 2 , derived 
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34A we derive slightly higher gas densities than for HH 34 
E. For HH 34 E we also have the two CS transitions, hence 
we attempted an LVG analysis by fixing the temperature 
to 10 and 15 K. As in the case of HH 7-11 there are no 
solutions for emitting sizes smaller than ~ 30" and the 
CS column density is derived assuming a volume density 
similar to the value derived from the HCO + . 

4.1.4. Summary of LVG analysis 

In summary, the LVG analysis leads to the following pre¬ 
liminary conclusions: (i) The gas densities of the clumps 
are always higher than those typical of a molecular cloud 
(~ 10 4 cm -3 ); (ii) HCO + (and CS for HH 7-11 and HH 34 
E) column densities derived from LVG are always larger 
than the beam averaged ones, consistent with the fact that 
the size of the clumps is smaller than the beam sizes, al¬ 
though in the case of the HH 7 11 clumps the picture is 
somewhat more complicated due to the close vicinity of 
the clumps; (in') The emitting sizes derived from the LVG 
analysis of CS seem to be larger than those derived from 
the HCO + analysis, indicating that the two species possi¬ 
bly arise from different gas components, a result already 
found by Girart et al. (2005) in their extensive study of 
the molecular gas ahead of HH 2. (iv) The temperatures 
are always well constrained and for all the clumps lower 
than 25 K, with HH7-11 being the warmest. Note that 
at low temperatures there is no clear pathway for the en¬ 
hancement of HCO + : we show in VW99 that for HCO + 
to reach column densities of > 10 14 cm -2 grain surface 
reactions and subsequent evaporation of the icy mantles 
must have occurred, as the main route of formation for 
HCO + is via the ion-neutral reaction of C + with H 2 O. 
If the temperature is < 100 K, mantle evaporation can 
only have occurred via photoevaporation. An alternative 
way might have been if the gas had been shocked in the 
past reaching higher temperatures and then cooled down. 
However during the time the gas cools back down to ~ 
10-20 K, HCO + should be destroyed again. Moreover, we 
show in VW99 that an enhanced UV radiation is in any 
case needed in order to ionize the gas enough to produce 
large abundances of C + . 

5. Discussion 

As explained in Section 2, our main criterion for extend¬ 
ing the survey to other molecular species was the strength 
of HCO + . Table [3] shows that most clumps that exhibit 
HCO + enhancement also show a rich chemistry, as pre¬ 
dicted by the VW99: methanol, when observed, seems to 
be at least one order of magnitude more abundant than 
what typically is found in dark molecular clouds or prestel- 
lar cores (see Table 5). Of the sulphur species, CS is ubiq¬ 
uitously detected but note however that CS is abundant in 
dark clouds and hence some of the emission may arise from 
the surrounding medium rather than the clumps, overall 
for those clumps with a size much smaller than our beam 
size. SO and H 2 S seem to be present in the majority of 


the clumps where HCO + is enhanced. Of the other sulphur 
bearing species, OCS, H 2 CS, and NS were only observed 
in HH 7 11 A. We will now analyse our results and look 
for a possible correlation between the chemical ‘richness’ 
of a clump, its vicinity to the HH shock and the strength 
of the radiation field impinging on the clump. 

In this section we attempt a simple analysis of the 
molecular species detected in our sample. In 5.1 we sum¬ 
marize the general physical and chemical properties of all 
the clumps. In 5.2 we discuss the effects of distance on 
the chemical processing ahead of HH objects by looking 
at the chemistry of different clumps ahead of the same ob¬ 
ject. However, it is important to note that the abundance 
of each molecular species will depend on i) the gas den¬ 
sity; ii) the radiation field ; Hi) the distance to the source 
of radiation; iv) the size of the clump, v) the age and, to a 
lesser extent, vi) the kinetic temperature of the gas. Some 
of these physical parameters have similar effects on the 
chemistry (e.g. the size and the gas density). Due to such 
degeneracy it is often difficult to disentangle the effects 
on the chemistry due to the different parameters. So, for 
example, for clumps ahead of the same object, where the 
radiation field is constant, different abundances may not 
be due solely to the different positioning of the clumps 
as it may be that the condensations differ in density, size 
and age. In 5.3 we briefly discuss the properties of the high 
velocity gas. 

5.1. Physical and chemical properties of clumps ahead 
of HH objects 

The physical and chemical properties of the gas ahead 
of the HH objects are different from those observed in 
prestellar cores or dark molecular clouds: they are denser 
than the latter but their small size and the high abun¬ 
dances of certain species, such as HCO + and methanol, 
make them unlikely candidates for pre-stellar cores. In 
Table 0 we compare the typical abundances of a pre¬ 
stellar core, L1544, with those derived for the HH 7 
11 A clump. The abundances of L1544 are taken from 
Jprgensen et al. (2002), while the HH 7 -11 A abundances 
were derived by dividing the averaged beam column densi¬ 
ties by the H 2 column density from the dust observations 
(IHatchell et al. 20051 by assuming a dust temperature of 
23 K IlChini et al. 200ijl . There are of course uncertainties 
associated with the derivation of the H 2 column densities 
(and hence of fractional abundances), such as i) the error 
bars associated with the integrated emission (see Table 10) 
and ii) the assumption of LTE which affects the molecu¬ 
lar abundances estimate. While the first one only leads 
to a 30% uncertainty, assuming LTE may lead to a larger 
uncertainty. The latter is still below the enhancement fac¬ 
tor found (i.e ~1 order of magnitude), although, for some 
species, the propagation of both errors may become rele¬ 
vant. Although the abundances ahead of HH 7 11 A are 
clearly enhanced with respect to those in L1544, the low 
abundance of CO is puzzling: while freeze out on the grains 
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is clearly responsible for its depletion in L1544, in the case 
of HH 7 11 A the UV radiation should have photoevap- 
orated all the ices so no significant amount of CO should 
be left on the dust unless it remains trapped in a mul¬ 
tilayered ice. An alternative possibility is that the con¬ 
densations we observe are transient and young: if the UV 
radiation photoevaporated the mantle while depletion was 
still occurring, there may have been not enough time to 
form CO. Viti et al. (2003) already observed that, in fact, 
to match the observations ahead of HH 2 the condensa¬ 
tions had to be young (< 1000 yrs) and in one of the best 
matching models, B10, CO has a gas phase abundance of 
2xl0 -5 at an Ay of 3.5 mags after evaporation. On the 
other hand it is also possible that the dust column density 
is overestimated: some of the material may not be associ¬ 
ated with the dense clumps and we may be overestimating 
the amount of dust responsible for the emission, perhaps 
because we have a too low an opacity. In short, it is hard 
to know what the absolute abundance of CO is and no 
strong conclusion should be drawn from Table El beside 
the fact that clearly some species are enhanced with re¬ 
spect to those in L1544. 

An enhancement of species such as HCO + , CH 3 OH 
and SO, in absence of shocks, can indeed be an indica¬ 
tion of active photochemistry. Jprgensen et al. (2002) do 
not give a fractional abundance for CH 3 OH for L1544 
but HCO + and SO are both less abundant by a factor 
of 5-10, although species such as CS have similar abun¬ 
dances. Interestingly DCO + also seems to be slightly en¬ 
hanced in HH 7-11 A with respect to the value found 
in L1544, although the DCO + /HCO + ratio is consis¬ 
tent with those expected in a dense, low temperature gas 
IRoberts et al. 2000l> . In fact a comparison between star¬ 
less cores and clumps ahead of HH objects may be mis¬ 
leading in that the anomaly of our clumps is that they 
are denser than a molecular cloud but much smaller than 
starless cores. They in fact show similarity with clumps 
detected by recent high resolution BIMA observations in 
L673 (IMorata et al. 20051 in that they are of similar sizes 
and they seem to have structures on scales of ~ 10 " or so 
but they have different chemistry because they are illumi¬ 
nated by the HH radiation. 

We are unable to compare clumps ahead of other re¬ 
gions with starless cores as we do not know their H 2 col¬ 
umn density but we compare the column densities across 
our clumps in Table Eland Figure 0 

Before attempting such comparison, it is important to 
note that, as discussed in Section 3 and as shown from 
our LVG analysis, we are uncertain about the size of all 
our clumps; hence for some objects, the beam size may be 
larger than the source size and for species such as HCO + 
(enhanced only in dense clumps subjected to the radi¬ 
ation) their averaged beam column densities may be a 
lower limit; however for species that are abundant also in 
the surrounding medium, such as CS and CO, their beam 
averaged column densities may be an overestimate of the 
emission from the clump. 


Table 5. Fractional abundances with respect to H 2 for 
HH 7-11 and L1544 


Molecule 

HH 7 
A 

-11 

X 

LI 544 

CO 

3.7(-5) 

5.5(-5) 

4.9(-6) 

CH 3 OH 

5.0(-9) 

9.8(-9) 

- 

HCO+ 

2.9(-9) 

3.8(-9) 

3.9(-10) 

H 2 CO 

l-3(-9) 

2.3(-9) 

- 

SO 

1.2(-9) 

2.0(-9) 

4.8(-10) 

CS 

8 . 6 (- 10 ) 

1.5(-9) 

8 . 6 (- 10 ) 

h 2 s 

T.l(-lO) 

l.l(-9) 

- 

HCN 

3.3(-10) 

3.8(-10) 

<4(-10) 

OCS 

3.0(-10) 

- 

- 

h 2 cs 

1 . 6 (- 10 ) 

- 

- 

DCO+ 

7.1(-11) 

8 . 6 (-ll) 

2T(-11) 


From Table El it is clear that all the clumps, except 
possibly JMG 99 G2, exhibit a similar chemistry: in or¬ 
der to eliminate the uncertainties due to the sizes, let us 
compare ratios rather than absolute column densities (see 
Figure 7); for example, HCO + /CH 30 H is remarkably con¬ 
stant within a factor of 2 for all objects, except for G 2 , 
where it is lower by one order of magnitude. Other ra¬ 
tios, such as HCO+/SO or SO/H 2 S also show remarkable 
consistency among the clumps although in G 2 the column 
densities of species characteristic of photochemistry are 
not as enhanced. In fact, chemically, it appears that while 
the clumps ahead of HH 7 11, HH 34 and HH 1 are of 
the same nature (confirmed also by the column densities 
derived from the LVG analysis), JMG 99 G1 and G2 have 
undergone different processing. 

Although we do not have an estimate for the temper¬ 
ature of JMG 99 G1 and G2, from our LVG analysis we 
note that HH 7 11 A clumps have a higher temperature 
than HH 34 and HH 1 clumps, therefore HH 34 and HH 1 
clumps seem to be the most quiescent, maybe because 
so far they only have been affected by the UV radiation. 
Indeed, the line widths of HH 34 are significantly smaller 
than those of HH 7-11 (Av(H 13 CO + ) ~ 0.4 and 1.2 km s _1 
for HH 34 and HH 7-11, respectively). HH 7 11 clumps, on 
the other hand, show sign of dynamical interaction with 
the HH objects, but such interaction seems not to have 
affected the chemical nature of the gas yet: the increase 
in temperature caused by the dynamical interaction is in 
fact still not sufficient to produce any shock chemistry and 
hence the chemical properties of HH 7-11 are very similar 
to those of HH 34 and HH 1 clumps. 

Note that, as also explained in Viti et al. (2004), the 
density enhancements we see ahead of HH 7-11 could not 
be produced by wind compression as the timescales would 
be too short for substantial freeze out to have occurred; 
as described in VW99 and Viti et al. (2003) depletion 
and subsequent hydrogenation on the grains must have oc¬ 
curred during the formation of the density enhancements, 
before the arrival of the HH object, in order for species 
such as HCO + and CH 3 OH to have formed. 
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It is tempting to draw an evolutionary scenario around 
these three different situations with HH 34 and HH 1 
clumps being the youngest and JMG 99 G1 and G2 being 
the oldest clumps. However, it is impossible to disentan¬ 
gle the effects of age, distance from the HH objects, and 
radiation, as discussed further in the next section. 

Despite the similar chemistry of most of the clumps 
and the uncertainties with the emission sizes, the (albeit 
small) differences in column densities (Table 3) may be 
significant. The three clumps ahead of HH 7-11 may be 
affected by several HH objects (Figure^, hence it is im¬ 
possible to determine which clump is the most irradiated; 
however, irradiation from several HH objects is consistent 
with the fact that HH 7 11 clumps show a richer chemistry 
than others. 

Although we do not here attempt any modelling of 
these clumps, it is interesting to compare the abundances 
of our clumps with the models computed in Viti et al. 
(2003). Note that these were detailed models of a partic¬ 
ular clump ahead of HH 2 hence we are not looking here 
for a perfect match. Viti et al. (2003) computed a large 
grid of models where both the radiation field and densities 
were varied. Although from the column densities in Table 
3 it seems that HCO + is higher ahead of HH 7 11 than 
ahead of HH 34 and NGC2264, from our LVG analysis it 
seems that the densities ahead of HH 7 11 are either lower 
or the same as those ahead of HH 34 and HH 1. Hence we 
may expect either that the radiation field from HH 7-11 
is higher, or that clumps ahead of HH 7 11 are younger. 
Figures 3 to 7 in Viti et al. (2003) show the column den¬ 
sities for HCO + , CH 3 OH , SO and H 2 CO as a function 
of time for an Ay ~ 3 mags. We compare the column 
densities in these figures with those derived from our ob¬ 
servations here for HH 7-11 and find that the best fitting 
models for HCO + , CH 3 OH and SO are B7 and B 8 which 
have a high radiation field of 1000 Habing and a density 
of 3xl0 5 cm -3 . The abundances of HCO + , CH 3 OH and 
SO may be less ahead of HH 34 by as much as 1 order 
of magnitude: although not shown in Viti et al. (2003) we 
compared two models from their grid where the density 
(3xl0 5 cm~ 3 ), age (< 100 yrs) and visual extinctions (~ 
3 mags) were the same but the radiation field varied from 
20 to 1000 Habing and we found that HCO + went from 
7xl0 12 to 5xl0 13 cm -2 , SO from 5xl0 14 to 2xl0 15 and 
CH 3 OH from lxlO 14 to 6xl0 14 cm -2 . However, at much 
later ages, a higher radiation field implies in fact a lower 
abundances of these species as the radiation field enhances 
but also destroys the gas: as explained in VW99 the effect 
of the irradiation is in fact to create a ‘wave’ of molecular 
formation where the abundances rise to enhanced values 
with respect to standard dark cloud chemistry but then 
fall again as photodissociation begins to dominate. It is 
difficult therefore to disentangle the effects of radiation, 
from those of age. 

Finally, it is interesting to note that in HH 7-11 A, 
the abundance of HCS + is at least one order of magnitude 
lower than that of H 2 CS. Again, from the lViti et al. 20031 
models we find that, as long as the density is greater than 


10 5 cm _3 then H 2 CS/HCS + is always > 10; but lower 
radiation fields (~ 20 Habing) yield to very high ratios 
(100-1000) while if % ~ 1000 Habing the ratio is ~ 10-50. 

5.2. Effects of distance on chemical processing ahead 
of HH objects 

In this section we attempt to single out the effect of the 
distance from other parameters under the assumption that 
the rich chemistry is mainly due to the UV impinging 
on the clumps. We therefore compare column densities 
among clumps belonging to the same region, where not 
only the radiation field from the HH shock is constant, 
but also where we assumed that the size of each clump 
is the same. In particular HH 34 may be the best region 
for this analysis as i) the four clumps show some chemical 
diversity and ii) unlike the HH 7 11 region there seems to 
be only one HH object impinging on the clumps. 

In the HH 34 region HH 34 C and E are the closest 
clumps to the HH object, with A and D 1.6 and 2 times 
further away (in the plane of the sky; the true distances 
from the HH object are of course unknown). From TableEl 
HH 34 C and E clearly show higher HCO + column den¬ 
sities than A and D, although the CO and CS column 
densities are similar, consistent with exposure to a higher 
UV field. Of C and E, E has the larger HCO + column den¬ 
sity, possibly because it has the higher derived gas density 
derived from the LVG analysis. Moreover, from Table 0 
SO and methanol from HH 34 E have higher column den¬ 
sities than the more distant HH 34 A. From this object, we 
also detect CsH 2 , not usually detected in quiescent dense 
clumps. This increase in chemical activity where the UV 
radiation field would be expected to be strongest is consis¬ 
tent with the UV chemistry, though a variation in density 
or age of the clumps cannot be ruled out as an alternative 
cause. 

5.3. The High Velocity component 

HH 7 11 A and JMG 99 G1 are the positions where the 
high velocity emission is brighter. We calculated the col¬ 
umn density ratio with CO for several species (from the 
13 CO spectra), i.e., the molecular fractional abundances 
with respect to the CO (see §4). Figure |H1 shows the 
fractional abundances with respect to the CO for several 
molecules as a function of the outflow velocity. 

From Figure 8 , it is clear that the HCO + relative abun¬ 
dance of the HH 7 11 high velocity component toward 
increases by more than an order of magnitude (a factor 
of ~ 30) between the slowest and fastest outflow compo¬ 
nents. This enhancement is also seen in H 2 CO and CH 3 OH 
for outflow velocities lower than ~ 6 km s~ x with respect 
to the systemic velocity. On the other hand, CS does not 
change its relative abundance. The HCO + abundance en¬ 
hancement at higher outflow velocities was first reported 
by Girart et al. (1999) towards the NGC 2071 molecular 
outflow. Recently, Dent et al. (2003) has found the same 
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Fig. 7. Ratios (in log 10 scale) of the column densities ol 
selected species for the low velocity gas. HH 7 B, HH 34 
C, HH 34 D and JMG 99 G1 are not shown due to lack of 
data. 



Fig. 8. Fractional abundance with respect to the CO of 
the high velocity gas for several molecules toward HH 7 
11 A. j/fiow = ULSR - «amb, where u a mb is the systemic 
velocity of the ambient gas. 


behaviour for the HCO + towards the high velocity gas 
around HH 2. However, HCO + is the only high density 
molecular tracer detected in the high velocity HH 2 com¬ 
ponent and this case could be a different scenario from 
NGC 2071 or HH 7—11, since in HH 2 the lack of other 
molecular high velocity gas other than HCO + can be ex¬ 
plained by the presence of strong UV field. In any case the 
clear enhancement of the HCO + abundance with respect 
to other molecular tracers with increasing high velocity in¬ 
dicates that it traces the molecular outflow well at higher 
velocities, whatever causes its enhancement. 

JMG 99 G1 not only has a spectral feature quite dif¬ 
ferent from the high velocity wings seen in HH 7-11 A but 
the chemical properties are also very different. In JMG 99 
G1 the CH 3 OH is by far the more abundant high density 
molecular tracer with a X[CH 30 H]/X[C 0 ]~ 2 - 4 x 10 -3 . 
This abundance is also found in strongly shock molecular 
clumps such those in LI 157 llBachiller et al. 200111 . and 
IRAS 21391+5802 UBeltran et al. 20021) . The HCO+on 
the other hand is almost one order of magnitude lower 
than in HH 7 11, although its abundance also increases 
with the outflow velocity. Even the spectral shape of the 
molecular emission in G1 is similar to the spectral shape 
of the lines detected in the L1157 clumps: Most of the lines 
shown are blueshifted (by 2-3 km s -1 with respect to the 
cloud velocity) and are significantly broader (Au ~ 4-5 
km s _1 ) than the lines from the other positions. Such emis¬ 
sion may in fact come from clumps similar to those seen 
along the L1157 and CB3 outflows IfBa,chiller et al . 2001 ’, 
IGodella et al. 199~9jl : Viti et al. (2004) found that such 
dumpiness may be a consequence of a pre-existing den¬ 
sity enhancement and of its interaction with the outflow. 


6. Conclusions 

We have performed a survey of molecular species ahead of 
several HH objects, where enhanced HCO + emission was 
previously observed. Recent modelling have supported the 
idea that such enhancement is due to UV radiation from 
the HH object lifting the ice mantles from the grains in 
the clumpy medium ahead of the jet, which then would 
drive a photochemistry. If this model is correct then other 
molecular species (such as CH 3 OH, SO etc) should also 
be enhanced. A recent molecular survey of the gas ahead 
of HH 2 seems to support such model (IGirart et al. 20051 . 
The aims of this present survey were to investigate: 1) 
whether the rich chemistry observed ahead of HH 2, and 
predicted by the VW99 model, was indeed characteristic 
of the gas ahead of all Herbig-Haro objects; 2) the chem¬ 
ical differentiation that may arise among clumps ahead 
of different objects, which may differ in radiation field 
strengths. 

We have selected 4 HH objects regions: in some of these 
regions enhanced HCO + emission had been detected from 
several chimps. Hence, in total we observed ten positions. 
For each position we searched for more than one species 
and, if possible, for more than one transition per species. 
We find that, chemically, the richest gas was that ahead of 
the HH 7 11 complex, in particular at the HH 7-11 A po¬ 
sition. In some regions we also detected a high velocity gas 
component. Beside estimating the averaged beam column 
densities for all the species observed, we have performed, 
where possible, LVG analyses to constrain the physical 
properties of the gas. The main conclusions from our LVG 
analysis are: 
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1. The gas densities of all the clumps observed are always 
higher than those typical of molecular clouds by at 
least a factor of 5, but probably higher. 

2. The derived gas temperatures are alway quite low (< 
15 K) apart for the clumps ahead of HH 7-11 where 
the maximum temperature does not in any case exceed 
25 K. 

3. HCO + , SO and CS column densities derived from our 
LVG and RADEX analyses are higher than the beam 
averaged ones. This indicates that the size of the emit¬ 
ting region is smaller than the beam size. 

4. The emitting size derived from the LVG analysis of CS 
are larger than those estimated from the HCO + anal¬ 
ysis. It may be that these two species arise from dif¬ 
ferent gas components, consistent with previous stud¬ 
ies f Girart et al. 2005 ). In addition, the HCO + abun¬ 
dance is enhanced with respect to the starless core in 
L1544, whereas the CS has similar abundances. All 
this suggests that at least a fraction of the CS emis¬ 
sion may arise from a cloud component or clumps not 
illuminated by he HH objects. 

By comparing column densities among the several 
clumps we also attempt a very simple chemical analysis in 
light of the photochemical model. We find that the abun¬ 
dances of the observed species ahead of the surveyed HH 
objects are higher than those found in a typical starless 
core and that by comparing the column densities among 
the different clumps we find that the clumps ahead of 
HH 7-11, HH 34 and HH 1 are of the same chemical na¬ 
ture, supporting the idea that the gas ahead of all HH 
objects exhibit a rich chemistry. 
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Table 6. Molecular line results for HH 7 -11 


Molecular 

Transition 


HH 7-11 A 



HH 7-11 X 



HH 7- 

-11 B 


Tmb 

(K) 

J TmbcIv 

(I< km s _1 ) 

ULSR 

(km s" 1 ) 

Aulsr 
(km s" 1 ) 

Tmb 

(K) 

J TmbcIv 
(K km s” 1 ) 

ULSR 

(km s" 1 ) 

Avlsr 
(km s _1 ) 

Tmb 

(K) 

J TmbcIv 
(K km s -1 ) 

vlsr 
(km s -1 ) 

Avlsr 
(km s -1 ) 

13 CO 2- 1° 

15.8(2) 

23.8(1) 

7.8(1) 

1.5(2) 

18.6(3) 

26.3(1) 

7.7(1) 

1.7(2) 

16.3(4) 

24.5(2) 

7.8(1) 

1.4(2) 

C ls O 2- 1° 

4.86(6) 

6.9(1) 

7.8(1) 

1.4(2) 

4.8(2) 

6.76(7) 

7.8(1) 

1.5(2) 





HCO+ 1-0° 

6.51(8) 

10.37(3) 

7.80(5) 

1.69(8) 

6.33(14) 

11.51(6) 

7.8(1) 

1.4(2) 

6.36(14) 

10.0(1) 

7.9(1) 

1.4(2) 

HCO+ 3-2 

9.5(3) 

18.08(18) 

7.87(1) 

1.80(2) 

8.5(6) 

17.7(4 ) 

7.92(2) 

1.94(5) 

10.0(6) 

18.1(4) 

7.8(1) 

1.6(2) 

H 13 CO+ 1-0 

1.91(3) 

2.43(2) 

7.92(1) 

1.20(1) 

1.71(5) 

2.21(3) 

7.89(1) 

1.21(2) 

1.92(6) 

2.44(4) 

7.97(1) 

1.20(2) 

H 13 CO+ 3-2 

1.0(2) 

1.08(8) 

8.03(4) 

1.03(9) 

0.71(18) 

0.95(9) 

8.11(6) 

1.26(12) 

1.1(3) 

1.3(2) 

8.21(8) 

1.2(2) 

DCO+ 2-1 

5.00(14) 

5.42(6) 

7.86(2) 

1.02(1) 

3.7(2) 

3.78(8) 

7.73(1) 

0.97(2) 





DCO+ 3-2 

3.5(2) 

3.45(8) 

7.84(1) 

0.91(2) 

2.1(3) 

2.48(15) 

7.75(3) 

1.10(7) 





HCN 3-2 

2.2(4) 

7.25(31) 

7.82(7) 

3.2(2) 

2.0(4) 

5.1(5) 

8.00(9) 

2.4(3) 





H 2 CO 2o,2-1o,i 

4.2(2) 

7.06(7) 

7.68(1) 

1.59(2) 

3.51(18) 

7.23(12) 

7.71(2) 

1.93(4) 





H 2 CO 3i, 3 —2i, 2 

4.2(1) 

6.42(3) 

7.86(1) 

1.42(1) 









CHaOH 2 0 -lo A+ 

2.62(5) 

3.60(3) 

7.83(1) 

1.29(1) 

2.38(7) 

3.53(4) 

7.84(1) 

1.40(1) 





CHaOH 2_i-l_i E 

1.93(5) 

2.67(3) 

fixed 

fixed 

1.74(7) 

2.60(4) 

fixed 

fixed 





CHaOH 2 0 -lo E 

0.36(5) 

0.49(3) 

fixed 

fixed 

0.29(7) 

0.44(4) 

fixed 

fixed 





CHaOH 2i-li E 

0.11(5) 

0.16(3) 

fixed 

fixed 

0.09(7) 

0.12(4) 

fixed 

fixed 





CHaOH 3o-2 0 A+ 

3.13(10) 

4.33(7) 

7.86(1) 

1.30(1) 

2.91(11) 

5.08(9) 

7.96(1) 

1.59(2) 





CHaOH 3-i-2_i E 

2.53(10) 

3.50(7) 

fixed 

fixed 

2.59(11) 

4.54(9) 

fixed 

fixed 





CHaOH 3o-2 0 E 

0.25(10) 

0.35(6) 

fixed 

fixed 

0.53(11) 

0.99(8) 

fixed 

fixed 





CHaOH 5o-4 0 A+ 

1.01(16) 

1.42(14) 

7.82(5) 

1.33(7) 

1.5(3) 

2.3(2) 

8.02(6) 

1.39(9) 





CHaOH 5-i-4_i E 

0.96(16) 

1.35(12) 

fixed 

fixed 

1.3(3) 

1.8(2) 

fixed 

fixed 





C 3 H 4 5o-4 0 

0.14(4) 

0.13(1) 

7.87(3) 

0.85(4) 









C 3 H 4 5i-4i 

0.14(4) 

0.13(1) 

fixed 

fixed 









C 3 H 4 5 2 -4 2 

0.08(4) 

0.07(1) 

fixed 

fixed 









CaH 2 2i j2 -1o,i 

1.43(12) 

1.33(3) 

8.10(1) 

0.86(2) 









CaH 2 3i, 2 -2 2 ,i 

0.53(10) 

0.69(8) 

7.87(1) 

1.21(10) 

0.50(11) 

0.65(8) 

7.96(1) 

1.21(11) 





CaH 2 41,4-30,3 

1.17(11) 

0.91(4) 

8.10(2) 

0.73(4) 









CS 3-2 

6.5(2) 

12.10(9) 

7.73(1) 

1.72(2) 

6.2(3) 

12.6(2) 

7.70(2) 

1.93(4) 

6.3(3) 

12.0(2) 

7.84(1) 

1.77(1) 

CS 5-4 

2.9(1) 

4.00(4) 

7.90(1) 

1.26(1) 









SO 3 2 -2i 

5.43(9) 

7.23(4) 

7.80(1) 

1.25(1) 

5.04(8) 

6.73(4) 

7.78(1) 

1.26(1) 





SO 65-54 

2.7(2) 

3.06(9) 

7.73(1) 

1.05(3) 

2.0(2) 

3.09(12) 

7.93(3) 

1.39(7) 





SO 7 6 -6 5 

1.1(3) 

1.7(2) 

7.61(8) 

1.4(2) 

0.8(3) 

1.1(2) 

7.8(1) 

1.2(4) 





H 2 S 1i,o-0 0 ,i 

2.8(2) 

4.17(11) 

7.73(2) 

1.41(4) 

2.1(3) 

3.9(2) 

7.77(5) 

1.75(11) 

2.2(3) 

3.5(2) 

7.92(5) 

1.50(9) 

H 2 CS 4 0 ,4-3o,a 

0.49(9) 

0.48(5) 

7.61(5) 

0.90(9) 









NS T/2-5/2 161.297 

0.16(7) 

0.18(3) 

7.94(7) 

1.05(8) 









NS T/2-5/2 161.298 

0.14(7) 

0.15(7) 

fixed 

fixed 









NS T/2-5/2 161.302 

0.14(7) 

0.16(3) 

fixed 

fixed 









OCS 7-6 

0.25(7) 

0.21(2) 

7.76(5) 

0.79(9) 









HCS+ 5-4 

0.21 

^ 0.10 












The 13 C0 TO, C ls O 1-0, HCO + 1-0 lines are clearly non-Gaussian (self absorption for the HCO + and strong high velocity emission for all of them). Line parameters 
derived from the zero (integrated intensity), first (line velocity) and second (line width) order moments of the emission in the 6.9-8.9 vlsr range (the velocity range of 
the H 13 CO + emission). 
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Table 7. Molecular line results for HH 34 E & A 


_ HH 34 A _ HH 34 E _ 

Molecular T M b J T M sdv wlsr Aw L sr Tmb J T M Bdv ulsr Au L sr 

Transition (K) (K km s -1 ) (km s _1 ) (km s _1 ) (K) (K km s _1 ) (km s _1 ) (km s _1 ) 


13 CO 2 1“ 

C ls O 2 1“ 

HCO+ 1-0 
HCO+ 3-2 
H 13 CO+ 1-0 
H 13 CO+ 3-2 
DCO+ 2-1 
DCO+ 3-2 
HCN 3-2 
H2 CO 2o,2-1o,i 
H2 CO 3i,3-2i,2 
CH 3OH 2 0 -l 0 A+ 
CH3OH 2_i-l_i E 
CH3OH 2 0 -l 0 E 
CH3OH 3o-2 0 A+ 
CH3OH 3_i-2_i E 
CS 3-2 
CS 5-4 
SO 3 2 -2i 
SO 6 s -5 4 
SO 7 6 -6 5 
C3H2 21,2-10,1 
C3H4 5o-4o 
H2CS 4o,4 _ 3o,3 
NS 7/2-5/2 161.297 


11.7(3) 

10.04(2) 

5.51(15) 

2.70(4) 

4.06(9) 

2.31(3) 

2.3(4) 

1.14(13) 

0.72(8) 

0.31(2) 

0.6(2) 

0.20(6) 

1.2(2) 

0.46(3) 

0.7(2) 

0.26(4) 

£ 0.75 

£ 0.20 

1.9(2) 

0.76(5) 

0.68(7) 

0.31(2) 

0.42(7) 

0.19(2) 

0.21 

£ 0.06 

£0.21 

£ 0.06 

£0.21 

£ 0.06 

2.5(3) 

1.37(6) 

1.60(7) 

0.80(2) 

£ 1.2 

£ 0.26 

£ 2.5 

£ 0.68 


8.60(9) 

0.65(9) 

8.70(1) 

0.46(1) 

8.49(1) 

0.54(1) 

8.51(3) 

0.46(5) 

8.63(1) 

0.40(3) 

8.3(1) 

0.3(1) 

8.54(1) 

0.37(3) 

8.55(4) 

0.37 6 

8.48(1) 

0.38(4) 

8.64(1) 

0.43(2) 

fixed 

fixed 


8.58(1) 

0.52(2) 

8.72(1) 

0.45(1) 


12.2(3) 

11.15(8) 

4.2(2) 

2.65(8) 

2.85(12) 

2.00(5) 

2.4(5) 

1.69(19) 

0.91(7) 

0.44(2) 

£ 1.0 

£ 0.22 

1.25(15) 

0.55(3) 

0.67(16) 

0.20(4) 

£ 0.75 

£ 0.20 

1.5(3) 

0.86(12) 

1.43(11) 

0.87(3) 

1.41(7) 

0.68(2) 

1.12(7) 

0.54(2) 

£0.21 

£ 0.06 

1.6(2) 

0.60(3) 

1.6(2) 

0.58(3) 

2.4(3) 

1.65(8) 

0.5(2) 

0.17(3) 

3.18(8) 

1.63(2) 

£ 1.2 

£ 0.26 

£2.5 

£ 0.55 

0.41(6) 

0.17(2) 

£ 0.19 

£ 0.05 

£ 0.36 

£ 0.09 

£ 0.45 

£0.11 


8.7(1) 

0.6(1) 

8.73(1) 

0.64(2) 

8.49(1) 

0.66(2) 

8.63(4) 

0.65(10) 

8.78(1) 

0.46(2) 

8.69(1) 

0.42(2) 

8.43(5) 

0.3(1) 

8.56(3) 

0.56(10) 

8.78(1) 

0.57(2) 

8.75(1) 

0.42(1) 

fixed 

fixed 

8.79(1) 

0.35(1) 

fixed 

fixed 

8.57(2) 

0.66(3) 

8.75(3) 

0.35(5) 

8.76(1) 

0.48(1) 

8.78(2) 

0.39(4) 


a The 13 CO 1-0 is non-Gaussian. Line parameters derived from the zero (integrated intensity), first (line velocity) and second 
(line width) order moments of the emission in the 8.30-9.20 vlsr range (the velocity range of C ls O emission). 

6 Weak line: fixed line width with the value of the DCO + 2-1 line. 


Table 8. Molecular line results for HH 34 C & D 


Molecular 

Transition 


HH 34 C 



HH 34 D 


Tmb 

(K) 

J TmbcIv 
(I< km s -1 ) 

ULSR 

(km s _1 ) 

Aulsr 
(km s -1 ) 

Tmb 

(K) 

J TmbcIv 
(K km s _1 ) 

ULSR 

(km s -1 ) 

Aulsr 
(km s" 1 ) 

13 co 2- r 

11.2(3) 

10.30(7) 

8.53(9) 

0.7(1) 

11.3(3) 

9.70(7) 

8.55(9) 

0.7(1) 

HCO+ 1-0 

3.79(9) 

2.95(2) 

8.51(1) 

0.73(1) 

2.78(9) 

2.22(3) 

8.57(1) 

0.75(1) 

HCO+ 3-2 

2.9(5) 

1.6(2) 

8.53(3) 

0.51(9) 

1.5(4) 

0.93(16) 

8.64(5) 

0.60(11) 

CS 3-2 

2.3(3) 

1.25(5) 

8.57(1) 

0.52(3) 

2.0(3) 

1.14(6) 

8.59(1) 

0.55(4) 


The 13 CO TO is non-Gaussian. Line parameters derived from the zero (integrated intensity), first (line velocity) and second 
(line width) order moments of the emission in the 8. 0-9.0 vlsr range (the velocity range of HCO + emission). 
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Table 9. Molecular line results for HH 1 


Molecular Tmb f TmbcIv wlsr Avlsr 

Transition (K) (K km s -1 ) (km s -1 ) (km s -1 ) 


13 CO 2 1“ 

C 18 0 2 1“ 

HCO+ 1-0 
HCO+ 3-2 
H 13 CO+ 1-0 
H 13 CO+ 3-2 
DCO+ 2-1 
DCO+ 3-2 
HCN 3-2 
H2CO 2o,2~1o,1 
CH 3 OH 2 0 -lo A+ 
CH 3 OH 2_i-l_i E 
CH 3 OH 2 0 -lo E 
CH 3 OH 3 0 -2 0 A+ 
CS 3-2 
SO 3 2 -2i 
SO 65-54 
SO 7 6 -6 5 

H 2 S li,o~ 0 o,i 


16.7(3) 

29.7(2) 

5.60(9) 

6.25(5) 

3.21(10) 

7.10(6) 

1.6(5) 

2.50(29) 

0.45(5) 

0.42(2) 

£ 0.33 

£ 0.15 

0.33(9) 

0.18(3) 

0.26(10) 

0.13(3) 

0.35(17) 

0.56(11) 

0 . 8 ( 2 ) 

1 . 2 ( 2 ) 

0.23(5) 

0.61(7) 

0.20(5) 

0.54(7) 

£ 0.18 

£ 0.13 

£ 0.21 

£ 0.11 

1.85(16) 

2.03(11) 

0.84(3) 

1.70(2) 

£ 0.6 

£ 0.28 

£ 2.0 

£ 0.93 

1.3(3) 

0.98(15) 


9.1(1) 

1 . 6 ( 2 ) 

9.1(1) 

1.3(1) 

9.21(1) 

2.07(2) 

9.12(8) 

1.45(20) 

9.31(2) 

0.87(5) 

9.19(4) 

0.57(13) 

9.18(7) 

0.45(13) 

9.14(16) 

1.5(3) 

9.18(9) 

1.5(3) 

9.20(10) 

2.4(2) 

fixed 

fixed 


9.09(3) 

1.03(8) 

9.05(1) 

1.92(3) 

9.01(6) 

0.73(14; 


a The 13 CO 1-0, C 18 0 1- lines are non-Gaussian. Line parameters derived from the zero (integrated intensity), first (line 
velocity) and second (line width) order moments of the emission in the 8.10-10.30 vlsr range (the velocity range of HCO + 
emission). 
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Table 10. Molecular line results for NGC2264G 


Molecular 

Transition 


JMG1 



JMG2 


Tmb 

(K) 

J Tyi&dv 
(K km s” 1 ) 

WLSR 
(km s -1 ) 

Avlsr 
(km s -1 ) 

Tmb 

(K) 

J Tmb(Iv 

(K km s _1 ) 

VLSR 
(km s _1 ) 

Aulsr 
(km s -1 ) 

13 CO 2- 1° 

1 . 0 ( 2 ) 

1.04(12) 

hxed 

hxed 

2.5(3) 

2.46(7) 

4.86(8) 

0.76(9) 


1.7(2) 

8.5(3) 

3.42(10) 

4.7(2) 

2.5(3) 

2.46(7) 

4.86(8) 

0.76(9) 


1.4(2) 

1.39(12) 

8.63(4) 

0.97(9) 

2.5(3) 

2.46(7) 

4.86(8) 

0.76(9) 

C ls O 2- 1 

0.61(18) 

0.65(9) 

4.83(7) 

1.01(18) 

1.04(19) 

0 . 86 ( 8 ) 

4.83(4) 

0.77(9) 

HCO+ 1-0“ 

0.7(1) 

2.30(6) 

1.6(4) 

3.0(5) 

0.4(1) 

0.33(3) 

4.7(2) 

0.9(2) 

HCO+ 3-2 6 

0.7(2) 

0 . 6 ( 2 ) 

4.9(1) 

0.9(3) 

2.8(4) 

2.77(14) 

4.6(1) 

0 . 8(1 ) 


0 . 8 ( 2 ) 

2.2(4) 

2 . 1 ( 2 ) 

2.5(6) 

2.8(4) 

2.77(14) 

4.6(1) 

0 . 8(1 ) 

H 13 CO+ l-0 d 

0.06(3) 

0.06(2) 

5.1(3) 

0.9(5) 

0.24(5) 

0 . 20 ( 2 ) 

4.64(4) 

0.78(7) 

H 13 CO+ 3-2 

£ 0.5 

£ 0.3 



£ 0.75 

£ 0.29 



DCO+ 2-1 

£ 0.15 

£ 0.07 



0.15(5) 

0.19(3) 

4.44(9) 

1 . 1 ( 2 ) 

DCO+ 3-2 

£ 0.18 

£ 0.08 



£0.27 

£ 0.08 



HCN 3-2 

0.69(14) 

3.87(18) 

2 . 10 ( 12 ) 

5.3(3) 

0.84(12) 

2.45(14) 

4.05(7) 

2.7(2) 

CH 3 OH 2 0 -l 0 A+ 

0.48(2) 

2.78(5) 

2.29(6) 

5.40(6) 

0.56(3) 

1 . 11 ( 2 ) 

4.43(2) 

1.83(3) 

CH 3 OH 2-i-l-i E 

0.35(2) 

2.03(5) 

hxed 

hxed 

0.39(3) 

0.75(2) 

hxed 

hxed 

CH 3 OH 2 0 -lo E 

0.09(2) 

0.53(5) 

hxed 

hxed 

0.07(3) 

0.13(2) 

hxed 

hxed 

CH 3 OH 2i-li E 

0.03(2) 

0.19(5) 

hxed 

hxed 

< 0.10 

< 0.03 



CH 3 OH 3o-2o A+ 

0.65(10) 

3.64(12) 

1.68(7) 

5.27(9) 

0.73(16) 

1.21(7) 

4.48(5) 

1.56(7) 

CH 3 OH 3_i-2_i E 

0.60(10) 

3.38(12) 

hxed 

hxed 

0.51(16) 

0.84(8) 

hxed 

hxed 

CH 3 OH 3 0 -2 0 E 

<0.28 

< 0.30 







CS 3-2 

0.91(14) 

1.10(14) 

4.88(5) 

1.15(11) 

1 . 6 ( 2 ) 

2 . 20 ( 12 ) 

4.58(3) 

1.28(2) 


0.78(14) 

4.5(3) 

1.88(15) 

5.4(3) 





SO 3 2 -2i 

0.35(5) 

2 . 21 ( 8 ) 

1.99(10) 

5.7(2) 

0.67(7) 

0.88(4) 

4.56(3) 

1.22(7) 

SO 6 5 -5 4 c 

0.15(7) 

0.80(10) 

1.6(3) 

4.7(5) 

0.36(8) 

0.57(6) 

4.22(7) 

1.4(2) 

SO 7 6 -6 5 c 

0 . 2 ( 1 ) 

0.72(18) 

1.6(3) 

3.5(5) 

0.3(1) 

0.17(5) 

4.3(1) 

0.5(2) 

H 2 S 1i,q-0o,i 

0 . 2 ( 1 ) 

0.9(2) 

2.8(4) 

3.9(9) 

0.3(1) 

0.5(2) 

4.3(2) 

1.5(4) 


0 The 13 CO 1-0 is fitted with three Gaussian components. One is forced to have the same velocity peak and line width as the 
C 18 0 emission. 

“ At the ambient velocity (~ 4.8 km s -1 ) the emission is missing. Computed the moments for the broader component using 
the CS velocity range (—0.6-4 .8 km s -1 ) 

6 Used two components fits as in CS 3-2. 
d Very marginal detection. 

c Marginal detection, so we used the -2.2 to 5.2 km s ' 1 velocity range to derive the integrated line intensities, velocity and 
line width. 



